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ABSTRACT 

ICR J00291-I-5934 is one of the seven accreting millisecond pulsars (AMPs) discovered so far. We 
report on the aperiodic timing and color analysis of its X-ray flux, using all the RXTE observations of 
the 2004 outburst. Flat-top noise and two harmonically related quasi-periodic oscillations, all of them 
at very low frequencies (0.01-0.1 Hz), were present in the power spectra during most of the outburst as 
well as a very high fractional variability (~50 %). These properties are atypical not only for AMPs but 
also for neutron star low-mass X-ray binaries (LMXBs) in general. There are instead some remarkable 
similarities with the variability observed in black hole systems, reinforcing the connections between 
these two types of LMXB, as well as some interesting differences. We note finally that the results of 
this paper are difficult to reconcile with interpretations where any break frequency of power density 
spectra scales inversely with the mass of the central object at an accuracy sufficient to distinguish 
between the masses of neutron stars and black holes in LMXBs. 

Subject headings: binaries: close — pulsars: individual (ICR J0029H-5934) — stars: neutron — 
X-rays: binaries — cataclysmic variables: individual (V709 Cas) 



I. INTRODUCTION 

Most low-mass X-ray binaries (LMXBs) contain 
rapidly spinning neutron stars (NSs), which only in a 
few cases become apparent as millisecond X-ray pul- 
sars. This growing sub-class of ac creting millisecond 
pulsars (AMPs; see IWijnandi l2005l for a review) con- 
sists of dim transients that typically go in outburst 
during a few weeks once every few years (so far the 
only exception is HETE J1900. 1-2455, still active at 
the moment of writing, more than a year after its dis- 
covery). AMPs are thought to be the evolutionary 
link between neutron star LMXBs and millisecond ra- 
dio pulsars, according to the so-called r ecycling scenario 
(|Bhattacharva fc van den Heuvell Il991l . for a review). 
Eight years after the first of the AMPs was discove red 
(SAX ,11808.4-3658. iWiinands fc van der Kli^ [1991 it 
is still unknown why these systems show pulsations and 
the rest of neutron star LMXBs do not, although some 
explanations have been put forward (e .g. screening of 
the magnet ic field. | Cumming et al.ll200il smearing of the 
pulsations. ^Titarc huk et aL 2002[ ). 

Aperiodic variability in the X-ray lightcurves of black 
hole and neutron star LMXBs (BH/NS-LMXBs), such as 
broadband noise or quasi-periodic oscillations (QPOs), 
offers a powerful tool to study the accretion flow in 
its inner region, where the bulk of the radiation is 
thought to originate and where gravitational fields are 
extreme (jron der Khsl I200a 120061) . Using this tool NS- 
LMXBs were classified according to their correlated 
timing and energy spectral behavior into Z and atoll 
sources (jHasinger fc: van der Klis'TgH). As noted sev- 
eral vears ago (e.g. Ivan d er Klis 1994), very similar phys- 
ical mechanisms are expected to drive the accretion pro- 
cess in BH- and NS-LMXBs, and the difficulties in find- 
ing unequivocal signatures of a black hole or a neu- 

^ Astronomical Institute "Anton Pannekoek" , University of Am- 
sterdam and Center for High-Energy Astrophysics, Kruislaan 403, 
NL-1098 SJ Amsterdam, Netherlands. 



tron star in the energy or power spectral properties of 
LMXBs are a natural consequence of this fact. Sev- 
eral attempts to find direct BH signatures have been 
made based on the s pectral prop erties at high energ ies 
(ISunvaev et al.lll991l : but see also ldi Salvo et al.|[2004l or 
iBarret et al.ll2000f ) or in the power spectral prop erties at 
high frequencies (ISunvaev fc RevnivtsevI [2000[ but see 
also iKlein-Wolt et al.ll2004f) . On the other hand, apart 
from pulsations and type I X-ray bursts, twin variable- 
frequency kilohertz quasi-periodic oscillatio ns (twin kHz 
QPO s) are now considered a NS signature (|van der KlisI 
|2006[ ). NS-LMXBs are also kno wn to be fainter than BH - 
LMXBs at radio wavelengths (iMig hari fc Fendeill200l l . 
and BH-LMXBs are g enerally fainter iii X-rays when 
they are in quiescence (jCarcia et al.ll200ll ). 

We present in this work the results of the aperiodic 
timing and X-ray color analysis of one of the seven 
known AMPs, ICR J00291-H5934. We show the ex- 
ceptional behavior of this source in its 2004 outburst, 
more similar with respect to aperiodic X-ray variability 
to a BH-LMXB than to the presently known NS-LMXBs: 
ICR J0029H-5934 was variable at very low frequencies 
with a very high fractional amplitude. However, previ- 
ously identified NS signatures were also present in the 
highest frequencies. In Section |2| we briefiy review some 
information on this source and in Sections |3| and |4| we 
describe our analysis and give its results, respectively. 
Section [5] contains our final remarks and conclusions. 

2. ICR J00291+5934 

On December 2'"^, 2004, ICR J00291-H5934 was dis- 
covered with IN TEGRAL dur ing one of the Galactic 
plane scans (|Eckert et al.ll20d^ . Follow up observations 
with RXTE revealed coh erent pulsations at 598.88 Hz 
(|Markwardt et all l2004bD , making it the fastest spin- 
ning AMP known. The pulse frequency exhibited a 
sinusoidal modulation with a period of 147.4 minutes 
(^2.45 hr) indicative of the orbital motion of t he pulsar 
in the binary system (|Markwardt et ai]|2004a| ). Optical 
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2004 iRupen et alJ l20dl and infrared (jsteeghs et al 



20041 ) counterparts were ide ntified and a recu rrence time 



of ~3 yr was suggested by iRemillardl ()2004D , based on 
RXTEs all-sky monitor dat a. Positive sp i n frequ ency 
derivatives were rep orted by iFalanga et alJ (|2005cf ) and 
iBurderi et al.l (|2005D . Several estimates of the distance to 
the source h ave been made, ranging frorn ^3 to ^10 kpc 
(iShaw et all [200 5: Gallow av et all 120051: iFalanga et all 
I2005cl: iJonkeret al.l 120051 : IBurderi et al.ll2005n . 

3. OBSERVATIONS AND DATA ANALYSIS 

We used all the pointed PCA observations of the 2004 
outburst of IGR J00291+5934, performed by RXTE be- 
tween December 3d. and 20th (see Table [Ij. An inter- 
mediate polar^ was also inside the PCA's field of view: 
V709 Gas (whic h has a spin period of ^313 seconds; 
Habcrl & Motch 1995 HNorton et al.lll999l : IFalanga et all 
2005b) . Hereinafter we refer to this system as the IP. 
All necessary caveats were taken in both color and tim- 
ing analysis in order to avoid contamination from this 
foreground source, as explained below. 

After filtering out unsuitable data according to the 
recommended criteria^ we extracted count rates from 
Standard 2 data, with 16 s time resolution, in the fol- 
lowing bands (in keV): A: 2.0-3.5, B: 3.5-6.0, C: 6.0- 
9.7, D: 9.7-16.0, and we computed soft (B/A) and hard 
(D/G) colors and intensity (A-fB-fG+D). The instru- 
mental background was subtracted according to the stan- 
dard faint-source background models and the colors and 
intensity normahzed to the Grab values nearest in time 
(jKuulkers et al.l[T994l) and in t he same PGA gain epoch 
(e.g.. Ivan Straaten et al.l [20 031. Using_WebPnnms and 
the spectrum reported by IFalang a et al.l (|2005H ). and 
taking into account the ~0.3 deg. offset between the 
IP and our source, we estimate a foreground 2-16 keV 
count rate from the IP of ~3.6 c/s/PCU (to be com- 
pared with the peak rate in our observations of ~76 
c/s/PCU ), which matches approximately the observed 
intensity during our last dataset (~5 c/s in PGU2), when 
IGR J00291-I-5934 had presumably ceased activity (see 
Section m Figure [1]). The effect of this contaminating 
flux is probably strongest on the colors, which are sensi- 
tive to subtle changes in the energy spectrum. For that 
reason we consider only the color-color diagram of the 
first half of the outburst, where the flux from the IP is 
less than 20% of the background subtracted total (see 
Section HT^ . 

GoodXenon and Event modes were used for the tim- 
ing analysis (see Table [iJ , rebinned in both cases to a 
1/8192-s time resolution (Nyquist frequency of 4096 Hz) 
and including only the counts in the ~2.5-30 keV range 
(absolute channels 5-71; this choice optimizes the S/N). 



^ Intermediate polars f Patte rson] 11994] . for a review) are a sub- 
class of accreting white dwarfs with a magnetic field strong enough 
to channel the accretion flow onto the magnetic poles (thereby 
emitting X-rays pulsed at the spin period) but unable to lock the 
rotation of the white dwarf with the orbital motion (as is the case 
for the so-called polars). Chance placed two semidetached binaries 
emitting X-ray pulses inside the PCA field of view, one of them 
containing a white dwarf and the other a neutron star. 

^ Elevation from the Earth greater than 10 degrees 
and pointing offset lower than 0.02 degrees; see PCA di- 
gest at http://hea sarc.nas a. gov/docs/xte/pca^news.html^ or, e.g., 
Ivan Straaten et al.l ll2U0cil) 



We performed fast Fourier transforms in 1024-s data 
segments, reaching thereby Fourier frequencies down to 
~0.001 Hz. Two observations were shorter than 1024 
s (90425-01-01-01 and 90425-01-01-03) so we performed 
and analyzed a 512-s FFT in each of them. No back- 
ground subtraction or dead-time correction were made 
prior to the FFTs and the Poisson noise power was sub- 
tracted &ran_yie resulting power density spectra, fol- 
lowing EinlmS (I2OOI : we first estimated the Pois- 
son noise using the expression proposed bv I Zhang et al.l 
([1995) and then (after inspecting the ^^2-4 kHz range 
and finding no unexpected features) shifted it to match 
the level between 2000-4096 Hz, where no intrinsic power 
is expected to be present, but only counting statistics 
noise (this shift was in all cases smaller than 0.06% of 
the previously estimated Poisson level). We then nor- 
mali zed the power spect ra in the so-called rms normaliza- 
tion ()van der Klislll995l ). For this normalization purpose, 
in order to take into account the foreground flux con- 
tributed by the IP, we used the count rates in the tail of 
the outburst (when the AMP reaches its quiescent state; 
set G; see Fig. [1]) to estimate the average background 
and foreground rate during the data segments used for 
timing. In order to improve the statistics we averaged 
observations successive in time and with power spectra 
consistent within errors, thereby dividing the outburst 
data into four sets (TablelU Fig[T]). The integrated power 
was calculated between 0.01 and 100 Hz in order to avoid 
contamination from the IP's variability (see Section [42| . 

We fitted the power spectra with a fit function com- 
posed of a sum of Lorentzians in the ^^Vmax repre- 
sentation'"*, excluding the first four frequency bins to 
preclude the IP's variability. In some cases, in or- 
der to avoid a meaningless negative coherence, Q was 
fixed to zero, which is equivalent to fitting a zero- 
centered Lorent zian. From now on, fol lowing previ- 
ous wo rk (Bcll oni et aT]l2002t W&n Straaten et all 120031: 
van der Klis.2006t lAltamirano et al.ll2005t iLinares et al.l 
we refer to these components as Li, where the L 
stands for Lorentzian and 'i' is the label identifying the 
component (such identification is given in Sections 14.2] 
and ]5.2p . Following this notation we call Li's character- 
istic frequency Vi and its coherence Qi. Five to seven 
Lorentzians were necessary to fit the power spectra, as 
described in detail in the next section. 

4. RESULTS 

4.1. Colors and intensity 

The time evolution of the count rate and colors can be 
seen in Fig.[TJ The lightcurve shows two clearly different 
slopes in a semi-logaritmic plot, indicating two different 
expo nential decay tim escales. As noted in previous works 
(e.g. IWijnai^ 120051 ) ■ this is common behavior among 
AMPs, with XTE J1807-294 the only kn own exception 
(|Linares et al.ll2005l : [FaTanga et al.ll2005af) . We fitted the 



* In this representation l]Belloni et al.l 120021 ^ . if ug is the 
Lorentzian's centroid frequency and A its HWHM (half width at 

half maximum), Vmax = + A'^ gives the characteristic fre- 

quency of the feature (near the centroid if it is narrow and near 
the half-width if it is wide). The quality factor Q = uo/2A is used 
in this representation as a measure of the coherence of the vari- 
ability feature. Its strength is given by the integral power (0— 00) 
whose square root, in the normalization we use, is the fractional 
rms amplitude of the variability. 
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first two parts of the (2.0-16.0 keV) lightcurve (A1+A2 
and B) and found e-folding times of 7.9 ± 0.3 and 2.45 ib 
0.05 da ys, similar to the ones reported in iFalanga et al.l 
()2005d ) from the INTEGRAL 20-100 keV lightcurve. 

Overall the source spectr um was hard compared 
to that of other AMPs (jvan Straaten et all (20051 : 
iLinares et al.l 120051) and it softened at both high (6.0- 
16.0 keV) and low (2.0-6.0 keV) energies during the first 
~9 days decay (see time evolution of soft and hard color 
in Fig.[T]). The resulting track in the color-color diagram 
in this first half of the outburst is shown in Figure [2l Af- 
ter the first 9 days (sets B and C in our analysis), the IP 
contributes with a significant fraction of the flux (>20% 
if all the flux in C comes from the IP) and the increase of 
the colors is likely due to that (intermediate polars are 
thought to be th e hardest X-ray emitt ing cataclysmic 
variables; see e.g. lde Martino et al.ll2004l ). 

4.2. Aperiodic timing 

The power spectra (Fig. [3]) on first inspection show an 
overall shape characteristic of low (luminosity) and hard 
(spectra) states of atoll sources, usually called island and 
extreme island states (where spectral hardness and total 
rms increase and luminosity and variability frequencies 
decrease when going from the island to the extreme is- 
land state). However, the overall (0.01-100 Hz) fractional 
rms amplitude is in the '^42-58 % range (Fig. [6l Ta- 
ble [2]), much higher than any value reported to date for a 
NS-LMXB (see e.g . lOlive et al.l[l99l iBarret et al.l[2000l: 
iBelloni et a"l]|2002f l. Five to six Lorentzian components 
were used to fit these power spectra. The best fit param- 
eters are displayed in Table [H The fit functions com- 
prise i) one broad Lorentzian at 0.03-0.05 Hz accounting 
for the flat-top noise and its break {Lb for "break"; see 
Sec. [3] for nomenclature) ii) two QPOs between 0.02 and 
0.05 Hz {Lqpo and its subharmonic, Lqpo/2) and iii) 
three zero-centered Lorentzians with characteristic fre- 
quencies covering the 0.5-0.8 Hz, 3-5.5 Hz and 25-70 Hz 
ranges (possibly Lh, Liow and L^, respectively; see Sec- 
tion [5?2] for further discussion). All components are sig- 
nificant at a > 3(7 level except Lqpo and Lqpq /2 in 
dataset B. 

The flat-top noise present in these observations can 
be alternatively described with a sum of two broad 
Lorentzians, one of them giving the first break and the 
other describing a further steepening above that break 
(a "wing" above the break in Frequency x Power repre- 
sentation; see Figure [5|). For completeness we report 
the parameters of those fits in Table [31 These fits are 
better from a statistical point of view (4.7cr, 2.3fT and 
3cr improvement in set Al, A2 and B, respectively, ac- 
cording to the F-test) and give a more detailed descrip- 
tion of the structure around the break, but are more 
difficult to compare with other sources (where usually 
single-Lorentzian or broken power-law breaks are mea- 
sured). Therefore, we use the frequencies of the first de- 
scribed six-Lorentzians fit function when comparing with 
frequencies of other sources (see Section [5]) . 

Assuming the IP did not drastically vary its behav- 
ior synchronously to the AMP's outburst phases, we can 
reject the possibility that an appreciable part of the vari- 
ability above 0.01 Hz comes from the IP. This is shown 
qualitatively by the lack of power above ^0.01 Hz in the 
power spectrum of the last dataset (Figure |4]), when the 



AMP is back to quiescence. Quantitatively, if e.g., Lqpo 
came from the IP it had an rms amplitude of ~27% of the 
IP's flux in set Al, and fell below the detection level in set 
C. A similar quantitative argument can be made for all 
variability components, as the integrated power between 
0.01 and 100 Hz in set C was zero within errors. We split 
the energy range into four contiguous bands (channels 5- 
9, 10-15, 16-25, 26-71: 2.5-4 keV, 4.5-6.5 keV, 6.9-10.6 
keV and 11.0-30.0 keV, respectively) and found no sig- 
nificant variations with energy in the strength, width or 
frequency of the detected QPOs or noise components. 
On the other hand a clear time evolution of both the 
QPOs and the broad-band noise is detected as shown in 
Fig. [51 the (0.01-100 Hz) fractional rms variability in- 
creased during the outburst (from 41.7 to 57.8%), while 
the frequencies of the QPOs and the other components 
decreased (from ~0.79 to ^0.50 Hz for Lh). 

5. DISCUSSION 
5.1. Exceptional island state 

Our analysis of the fastest spinning AMP, 
IGR J00291-f 5934, reveals the exceptional behav- 
ior of this source in its 2004 outburst compared to 
the rest of NS-LMXBs: very strong X-ray variability 
at very low Fourier frequencies. Its power spectra 
showed flat-topped noise (i.e. approximately flat below 
a certain "break" frequency and decreasing above that) 
with extremely low break frequencies (^0.04 Hz), as 
well as very high integrated fractional rms variability 
(~50 %). These break frequencies are about two orders 
of magnitude lower than those observed in Z sources 
and about one order of magnitude lower than those 
previously found in extreme island (low-hard) states 
of at oll sources (e.g, IBarret et al.l 120001 : IBelloni et al.l 
|2002|) ^. The measured fractional rms amplitude of the 
variability is the highest found so far for in a NS-LMXB. 
Both the break frequencies and fractional variability 
observed in IGR J00291-I-5934 are instead similar 
to those detected in low-hard states of BH-LMXBs 
(though the frequencies are low even for the lowest 
BH-LMXB range; see next section for a more detailed 
comparison). It h as been noted (jvan Straaten et all 
120051: ILinares et al] I2005D that AMPs hardly "make 
it to the soft state" and instead remain in (low-hard) 
island or extreme island states. The (2.5-25.0 keV) 
peak luminosity of IGR J00291-I-5934 (between 0.5 and 
5.5 % of the Ed dington luminosity according to the 
flux measured bv lGallowav et al]|2005l and the different 
distance estimates, see Sec. [3]) was similar to that of 
other transient NS-LMXBs when they show "canonical" 
extreme island states, which seems to indicate that 
the exceptional state of our source was not purely a 
consequence of an exceptionally low accretion rate. 

The lack of high frequency variability (>100 Hz) dur- 
ing the outburst of IGR J00291+5934 might be related to 
a combination of a dynamically important magnetic field 
(needed in AMPs in order to produce the observed pulsa- 
tions) and a fast spin (spin frequency of ~600 Hz in the 
case of IGR J00291+5934) that could prevent the accre- 
tion disk from reaching the innermost regions, where the 
highest dynamical frequencies are present. This could in 

5 A reanalysis of GS 1826-238 data in IBarret et aP |[2000D led 
to a break frequency above 0.1 Hz 
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turn be related to the strong overall va riability that we 
measured. If, as has been argued (e.g. IChurazov et al.l 
I2OOII: iGilfanov et"al] l2003l ). the disk sets the variabil- 
ity frequencies and the comptonizing medium (corona, 
boundary layer, spreading layer, ADAF or similar) sets 
the variability amplitude, then a disk truncated at a 
large radius could provide both low dynamical frequen- 
cies and strong fractional variability, as its flux contri- 
bution would be relatively low and hence so would be 
the decrease that this flux would cause in the relative 
amplitude of the variabilit y. We note, however, that in 
standard descriptions (e.g. lGhosh fc Lamblll978f l the in- 
ner disk radius does not depend on the spin rate, and if 
the so-called propeller mechanism is responsible for clear- 
ing out the inner disk annulus this woul d be in contra- 
diction with the rep orted spin- up (Fala nga et al.l[2"005cl : 
iBurderi et al.l 120051 ) . as the expelled matter would re- 
move angular momentum from the neutron star. The 
interaction between a rapidly spinning magnetized neu- 
tron star and an accretion disk is a complex problem. 
Most of the work so far concerns the "classical" slo w X- 
ray pulsars (jLamb et al.lll973l : iGhosh &: Lambl[T978D and 
only recently the AMP case has received more attent ion 
(jPsaltis fc Chakrabartvlll999l : iRappaport et al.ll2004[ ). 

An anticorrelation in the flat-topped noise of Cyg X- 
1 was discovered by Belloni & Hasingcr (1990) linking 
the rms level of the flat top with the break frequency: 
the lower the break frequency, the higher is the flat-top 
(in rms normalization). In Figure [7] (top) we see how 
IGR J00291-I-5934 extends this relation to the lowest fre- 
quencies and the highest rms levels, even when compar- 
ing with the BH system Cyg X-1 fsee lBelloni et al.ll2002i 
for an update of this correlation and comparison with 
other neutron star systems, which are all at frequencies 
>0.1 Hz, and the BH-LMXB XTE J1118-I-480, which 
is contiguous with but still at slightly higher frequecies 
than our NS-LMXB IGR J00291+5934). 

5.2. IGR J00291+5934 vs. BH-LMXBs 

The similarities between IGR J00291-H5934 and BH- 
LMXBs in low-hard states are illustrated in Figure [3 
where a power spectrum of IGR J00291-I-5934 is com- 
pared to one of XTE J1118-H480 and one of XTE JI550- 
564 (e.g. lRevnivtsev et al.ll2000l : ICui et al.lll999l respec- 
tively). The flat-top noise has very similar power, break 
frequency and QPOs. There are also obvious yet in- 
teresting differences: as previously noted in a study 
of low states of BH/NS-LMXBs the NS power spectra 
show mor e power at high frequencies ( >10 Hz) than the 
BH ones (jSunvaev fc Revnivtsev|[2000D . In this extreme 
case, however, the differences between the power spec- 
tra of IGR J00291-h5934 and XTE J1118+480 start at 
~0.2 Hz, where the NS power spectrum rises and the 
BH one starts to fall (see Fig. [8]). The power spectrum 
of XTE J1550-564 has a very similar shape to that of 
IGR J00291-I-5934 and XTE Jl 118-1-480, but shifted to- 
wards higher frequency and lower rms. It also shows a 
prominent bump at ^I-IO Hz (similar to what is seen in 
IGR J0029H-5934 at slightly lower frequency) followed 
by a steep d ecline above ^10 Hz (ch aracteristic of BH 
systems; see iKlein-Wolt fc van der K lis 2006 for an ex- 
tensive study and furt her discussion). I n contrast to the 
ones studied in ,Sunvaev fc RevnivtsevI (2000.) this NS- 
LMXB shows no significant power above ~100 Hz. 



We found two low-frequency QPOs at ~0.02 and 
~0.04 Hz, similar to the ones seen in BH- and other 
NS-LMXBs in low- luminosity states (jCasella et al.ll2005l: 
Ivan Straaten et all [20051. Their centroid frequency ra- 
tios were 2.06±0.06 and 1.91±0.06 in Al and A2 re- 
spectively, indicating an harmonic relation. In general 
this kind of QPO is associated with low-luminosity states 
and therefore it requires a relatively low mass accretion 
rate in the inner regions of the disk. Several attemps 
have been made to investigate the nature of these low- 
frequency QPOs, e. g. in the context of shot-noise mod- 
els jVi khlinin et al. 19941 ) or thermal- viscous instabilities 
(|Chen fc Taam. 19941 ). In any case our work shows clearly 
that neither the presence of a solid surface nor a magnetic 
field affects this phenomenon so that the same physical 
mechanism must be at work in NS and BH systems in 
order to produce this common feature: harmonically re- 
lated QPOs superposed on a break in the flat-top noise 
a t very low frequencies. 

IWijnands fc van der Klid (jl999 l) found a correlation 
between the break frequency and the frequency of the 
"hump" present in the power spectra of atoll, Z sources, 
black hole candidates and AMPs {Lh or sometimes an as- 
sociated QPO, Llf)- We find that the frequency of the 
first zero-centered Lorentzian detected above the break 
in the power spectra of IGR J00291-I-5934 also satisfies 
this WK correlation (Figure [71 bottom) which leads us 
to identify it as Lh- The fact that vqpo and vqpo/2 
fall much below this relation (on top of the break as dis- 
cussed) suggests that they are not directly related to the 
low frequency QPOs {L^p, Li^p/2) measured sometimes 
on top of the "hump" . Interestingly, the same happens 
for XTE J1118-h480 and XTE J1550-564: the QPOs on 
top of the break fall clearly below the WK correlation, 
whereas Lh follows it. 

The relations between the characteristic frequencies 
of several variability components of Z and atoll sources 
and the frequency of the upp er kilohertz QPO, y ,, , were 
studied and synthesized by Ivan Straaten et all ()2003D 
(Figure |9l). Shifts in these frequency-frequ ency corre- 
lation s were found in two AMPs by van Stra aten et al.l 
l|2005l ) and iLinares et all (|2005l in SAX J1808.4-3658 
and XTE J1807-294, respectively), in a way that was 
best explained by the upper and lower kHz QPOs having 
frequencies {vu, i^i) lower than the ones measured in atoll 
sources by a factor '^1.5. The component Lfow plotted in 
Figure [9] is the broad component present below L„ (the 
component with the highest characteristic frequency) in 
extreme island states of low-lu minosity atoll sources (in- 
cluding bursters and AMPs; Ivan Straaten et al] |2002| . 
2003, 2005,). We tentatively identify the two highest 
frequency variability components in IGR J0029H-5934 
as Lu and Lioyj, as these are the only broad variable- 
frequency components detected in extreme island states 
of NS-LMXBs above Lh, and plot the frequencies of 
Lb, Lh and Lgow against that of L„ (Figure [9K). As 
can clearly be seen in Figure [51 the frequencies of 
IGR J00291-I-5934 also vary in correlation to one another, 
but with slopes much less steep than those found in atoll 
sources and AMPs at higher frequencies. In order to 
compare IGR J00291-I-5934 with the two shifted AMPs, 
wc multiply and veow by a trial factor of 1.5 (Fig- 
ure Eb). The tracks traced by IGR J00291+5934 seem 
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to connect better with the ones at higher frequencies in 
the shifted version of the frequency- frequency plots, but 
of course the difference in slopes remains. We conclude 
that a radical change in the slope of these correlations 
appears to occur at frequencies i>u ^100 Hz. It is unclear 
at this point to what extent this is universal to LMXBs 
or specific to IGR J00291+5934. It will be interesting 
to study the behavior of our source in future outbursts 
to see if at I'u > lOOHz its components follow the usual 
correlations or not. 

5.3. Mass scaling for break frequencies 

This work also sheds some light on another interest- 
ing issue regarding broad-band noise in power spectra. 
Break frequencies are often used to estimate masses of 
accreting black holes with the argument that timescales 
in the accretion disk scale linearly with the mass of 
the central object and theref ore the observ ed frequencies 
should scale inversely (e.g., lUttlev et al.. 2002). From 
the theoretical point of view this argument is weak- 
ened by our poor knowledge of what produces the 
broad-band noise and which timescale sets its break 
frequency (viscous, dynamical, relate d w ith the size of 
the disk, with shot lifetime; see e.g. iPri edhorsky et all 
1979; Manmoto et al. 1996; Poutanen fc Fabiaiirfl999l : 
Merloni et al.i.200O: .Churazov et al.,.2001i) . From the ob- 
servational point of view the application of this argu- 
ment using break frequencies {vb) of the flat-top noise 
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(Lb) observed in low-hard states of LMXBs is clearly 
challenged by IGR J00291-I-5934, as it harbors a neu- 
tron star and shows break frequencies lower than the 
ones observed in most of the more massive stellar-mass 
black holes. Although vi, of Cyg X-1 in the low-hard 
state was used in the past to infer masses of super- 
massive black holes, recent advances in active galactic 
nuclei (AGN) X-ray timing indicate that most AGN are 
in a state similar to the high-soft state of BH-LMXBs 
(jUttlev fc McHardvll2005f ). Any power spectral feature 
that is used in such mass-scaling arguments should i) 
have a rather stable frequency (like perhaps the hecto- 
hertz L orentzian in NS systems.lvan Straaten et "ani2002t 
see also lPottschmidt et al.ll2003D or ii) be accuratelv cor- 
rected for drifts due to some parameter varyi ng in and be- 
tween systems (e.g. mass accret ion rate; iMcHardv et al.l 
[200l lUttlev fc McHardvl 120051) or a combination of i) 
and ii). In any case the power spectra observed in 
IGR J0029H-5934 refute a simple mass-frequency rela- 
tion for most of the variability components present in 
LMXBs. 
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TABLE 1 
Log of the observations. 



Set 


ObsID'^ 


Date' 


Detectors^ 


Count Rate ^ 
(c/s) 


1024-s PDS^ 


Total / Bkg. ^ 
(c/a) 


Al 


90052-03-01-00 


0.9 


3.0 


308.9±0.0 


1 


232.9 / 61.2 




90052-03-01-14 


2.7 


3.0 


262.5±15.0 


2 






90052-03-01-04 


3.6 


3.0 


240.4±4.7 


12 






90052-03-01-05 


4.0 


3.2 


256.4±33.7 


8 






90052-03-01-06 


4.6 


4.0 


300.4±8.9 


3 






90425-01-01-08 


4.9 


3.3 


236.5±29.0 


7 






90425-01-01-02 


5.1 


5.0 


347.S±6.3 


3 






90425-01-01-01 


5.2 


5.0 


351.1±0.0 


_ 






90425-01-01-03 


5.3 


5.0 


369.0±0.0 


_ 






90425-01-01-09 


5.6 


3.0 


209.1±4.4 


3 






90425-01-01-07 


5.9 


3.5 


225.5±31.6 


12 






90425-01-01-10 


6.5 


2.9 


180.5±53.6 


14 






90425-01-01-12 


6.7 


3.0 


188.2±11.5 


2 




A2 


90425-01-01-110 


7.0 


3.3 


192.6±41.2 


13 


166.6 / 59.2 




90425-01-01-11 


7.2 


3.0 


166.9±1.1 


2 






90425-01-01-15 


7.3 


3.0 


168.3±2.1 


2 






90425-01-01-13 


7.5 


3.2 


174.3±21.3 


15 






90425-01-02-07 


7.6 


3.0 


156.7±0.0 


1 






90425-01-02-08 


7.7 


3.0 


174.6±15.2 


2 






90425-01-02-000 


7.9 


3.2 


167.3±22.2 


12 






90425-01-02-00 


8.1 


3.0 


148.2±2.4 


4 






90425-01-02-15 


8.2 


3.0 


150.3±2.4 


2 






90425-01-02-04 


8.3 


3.0 


144.6±3.8 


2 






90425-01-02-05 


8.4 


3.0 


142.6±2.6 


3 






90425-01-02-16 


8.5 


3.0 


141.5±4.5 


10 




B 


90425-01-02-02 


8.7 


3.0 


127.9±0.0 


1 


94.1 / 53.9 




90425-01-02-030 


8.9 


3.0 


126.6±5.4 


12 






90425-01-02-03 


9.1 


3.0 


117.5±3.7 


5 






90425-01-02-09 


9.2 


3.0 


112.6±3.5 


2 






90425-01-02-100 


9.4 


2.8 


102.3±14.0 


15 






90425-01-02-10 


9.6 


3.0 


126.0±0.0 


1 






90425-01-02-11 


9.7 


3.0 


99.9±0.0 


1 






90425-01-02-120 


9.9 


2.8 


91.8±9.2 


12 






90425-01-02-12 


10.1 


2.8 


81.9±10.7 


5 






90425-01-02-13 


10.3 


3.0 


85.9±2.7 


3 






90425-01-02-14 


10.4 


2.6 


77.7±13.9 


11 






90425-01-02-17 


10.6 


2.0 


65.7±7.2 


2 






90425-01-02-06 


11.5 


3.0 


72.5±3.S 


5 






90425-01-02-01 


11.6 


2.0 


52.1±0.0 


1 






90425-01-02-18 


11.9 


3.3 


74.0±17.6 


10 






90425-01-02-20 


13.2 


3.0 


58.1±0.0 


1 






90425-01-02-25 


14.1 


4.0 


72.4±0.4 


2 


63.6 / 63.6 




90425-01-02-26 


14.4 


2.0 


40.0±0.0 


1 






90425-01-02-24 


14.5 


4.0 


77.4 ±3.7 


3 






90425-01-02-19 


14.5 


3.0 


72.8±0.0 


1 






90425-01-02-27 


14.6 


5.0 


91.5±0.0 


1 






90425-01-03-00 


14.9 


3.3 


63.5±8.1 


9 






90425-01-03-01 


15.3 


3.0 


56.1±2.5 


8 






90425-01-03-02 


16.0 


4.0 


72.7±1.7 


2 






90425-01-03-03 


16.2 


3.0 


55.2±1.5 


6 






90425-01-03-04 


17.3 


3.4 


66.0±10.1 


7 






90425-01-03-05 


17.9 


4.0 


71.9±2.1 


6 






90425-01-03-06 


18.2 


3.0 


54.8±2.0 







Observation start date in days since MJD 53341 (December 2"'', 2004, the discovery date) 
'^Average number of active detectors 

''Average and standard deviation of the ~2.5-30 keV count rate, including all active detectors and not corrected for background 
''Number of power density spectra extracted from the observation, each of them 1024 s long 

^Set averages of total and background ~2.5-30 keV count rates (background estimated from set C; see Section[3ll 
*Set C contains mainly background/foreground observations. 
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TABLE 2 

Best fit parameters from 6-Lorentzian model and integrated rms variability 



Set 


Parameter 


Lt 




Lqpo 


Lh 




Lu 


xVd.o.i. 


rms (%) 


Al 


I^max(Hz) 


(4.9±0.2)xl0-2 


(2.15±0.05)xl0-2 


(4.43±0.10)xl0-2 


0.79±0.03 


5.5±0.2 


70±10 


357/325 


41.7±0.4 




rms (%) 


24.6±0.3 


4.8±0.8 


4.5±0.7 


23.4±0.2 


22.1±0.2 


15.5±0.4 








Q 


(4±4)xl0-2 


5.5±3.2 


6.2±2.2 


(fixed) 


(fixed) 


(fixed) 






A2 


I^max{Hz) 


(3.9±0.1)xl0-2 


(2.11±0.09)xl0-2 


(4.03±0.06)xl0-2 


0.66±0.02 


4.4±0.3 


37.6±7.2 


335/326 


47.3±0.6 




rms (%) 


29.3±0.3 


5.5±1.0 


4.4±0.6 


27.1±0.3 


23.5±0.4 


15.6 ± 0.7 








Q 


(fixed) 


5.1±2.5 


16±6 


(fixed) 


(fixed) 


(fixed) 






B 




(2.8±0.1)xl0-^ 






0.50±0.03 


2.7±0.4 


25.0±7.7 


418/332 


57.8±1.2 




rms (%) 


37.5±0.3 






34.4±0.4 


25.9±1.0 


21.2±1.4 








Q 


(fixed) 






(fixed) 


(fixed) 


(fixed) 







TABLE 3 

Best fit parameters from 7-Lorentzian model 



Set 


Parameter 


Lib 


Lqpo/2 


Lqpo 


L2b 


Lh 


^£ow 


Lu 


X^/d.o.f. 


Al 




(1.6±0.3)xl0-^ 


(2.25±0.05)xl0-^ 


(4.36±0.08)xl0-^ 


(7.8±0.8)xl0~^ 


0.71±0.03 


5.3±0.2 


66.3±9.3 


327/322 




rms (%) 


12.8±2.1 


6.5±2.1 (2.3(7) 


8.3±1.6 


17.8±1.7 


24.0±0.2 


22.3±0.2 


15.6±0.4 






Q 


0.53±0.17 


3.6±1.6 


2.9±1.0 


0.44±0.11 


(fixed) 


(fixed) 


(fixed) 




A2 


I^max(Hz) 


(1.8±0.4)xl0-2 


(2.19±0.08)xl0-2 


(3.98±0.05)xl0-2 


(6.5±0.9)xl0-2 


0.63±0.03 


4.3±0.3 


36.7±6.9 


323/322 




rms (%) 


18.7±3.4 


5.8±2.0 (1.8cr) 


6.1±0.9 


21.0±3.1 


27.6±0.3 


23.7±0.4 


15.7±0.7 






Q 


0.3±0.1 


5+^ 


7.8±2.5 


0.3±0.1 


(fixed) 


(fixed) 


(fixed) 




B 


I^maxiHz) 


(1.8±0.3)xl0-2 






(5.0±1.3)xl0-2 


0.56±0.03 


3.1±0.5 


27±^^ 


403/330 




rms (%) 


29.5±4.1 






25.2±4.7 


34.5±0.7 


24.9±1.1 


20.7±1.5 






Q 


(fixed) 






(fixed) 


(fixed) 


(fixed) 


(fixed) 





Within this function two Lorentzians are used to fit the break in the fiat-top noise, Lu, and L26 (see Sec. 14.21 1 
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Fig. 1. — Time evolution of colors and intensity during IGR J00291+5934's December 2004 outburst. The upper panel shows the 
lightcurve from both ASM (1.5-12 keV, daily average; open circles) and pointed PCA (2-16 keV, observation average; filled symbols) data. 
Vertical dotted lines show the borders between averaged data sets Al, A2, B and C. The increase in colors after day ~10 is likely due to 
foreground contamination (see text). 
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Fig. 2. — Color-color diagram of data sets Al and A2, where the colors are not yet dominated by the intermediate polar in the field of 
view. The upper right point corresponds to the first pointed observation and the line connects observations successive in time. 
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Fig. 3. — Power spectra in Power x Frequency vs. Frequency representation. Also shown are the respective fit functions and their 
Lorentzian components. The spin frequency of the neutron star shows as a spike at ~600Hz. 



12 




0,01 0.1 1 10 100 1000 10 

Frequency (Hz) 



Fig. 4. — Power spectra of the background/foreground dominated set, in Power X Frequency vs. Frequency representation. The spin 
frequency of the white dwarf in the field of view is visible at ~0.003Hz (312. 8s), as well as its second harmonic. Both periodicities have 
been already detected in this system in the 0.1-2.4 keV i?05/lT energy band l IHaberl fc Motchlll99"5l : [Norton et al.|[T999f) 
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Fig. 5. — Zoom into the low-frequency part of the power spectrum of set Al, showing an alternative function which uses two broad 
lorentzians (dashed lines) to fit the break. See Section |3] for details. 
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Fig. 6. — Integrated (0.01-100 Hz) fractional rms amplitude and characteristic frequency of the "hump" in the power spectra {L}^) for 
the three first datasets. Error bars in the X axis indicate the time interval of each dataset. 
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Fig . 7. — (Top:) Correlation between rms level at the break of the power spectra and the frequency of this break (afterlBelloni &: Hasingeil 
119901 . BH90: Cyg X-1, grey squares). (Bottom:) Correlation between the frequencies of Lf, and Lj^, after [Wiinands & van 
The lines indicate the best power law fits to the original correlations and the points that come from our work are shown in black. 
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Fig. 8. — Power spect ra of the accreting mill isecond pulsar IG R J00291+5934 and of two black hole low-mass X-ray binaries in the low 
state (XTE J1118-I-480. IRevnivtsev et al.il2000l : XTE J1550-564. ICui et al.|[T999t) . The power spectrum of IGR J00291+5934 corresponds 
to sets Al, A2 and B of our analysis. The energy band of XTE J1550-564's power spectrum (~2-13 keV) is different than the one used in 
the other two cases (~2.5— 30 keV). The fractional rms variability (0.01-100 Hz) is shown on each plot. T he m imetism between the neutron 
star and the black hole candidates is obvious, although there are also interesting differences (see Section |5. 2^ . 
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Fig. 9. — Frequency-frequency correlations in neutron star low-mass X-ray binaries. IGR J00291-I-5934 (black symbols at i/u ^100 Hz) 
is compared with two accreting millisecond pulsars showing shifts in these correlations (SAX J1808.4— 3658 and XTE J1807-294, open 
symbols) and with non-pulsing atoll sources (grey symbols, after van Straaten ct al. 2002). In a) the characteristic frequencies of power 
spectral components (identified by the symbols on the left) are plotted versus the frequency of the upper kilohertz QPO (or Lu). In b) Uu 
and ui of SAX J1808.4— 3658 and XTE J1807-294 have been multiplied by the reported shift factors, whereas a factor of 1.5 has been used 
for IGR J00291+5934. 



